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Since its discovery in 2004, graphene, a two-dimensional hexagonal carbon allotrope, has gen- 
erated great interest and spurred research activity from materials science to particle physics and 
vice versa. In particular, graphene has been found to exhibit outstanding electronic and mechan- 
ical properties, as well as an unusual low-energy spectrum of Dirac quasiparticles giving rise to 
a fractional quantum Hall effect when freely suspended and immersed in a magnetic field. One 
of the most intriguing puzzles of graphene involves the low-temperature conductivity at zero den- 
sity, a central issue in the design of graphene-based nanoelectronic components. While suspended 
graphene experiments have shown a trend reminiscent of semiconductors, with rising resistivity at 
low temperatures, most theories predict a constant or even decreasing resistivity. However, lattice 
field theory calculations have revealed that suspended graphene is at or near the critical coupling 
for excitonic gap formation due to strong Coulomb interactions, which suggests a simple and 
straightforward explanation for the experimental data. In this contribution we review the current 
status of the field with emphasis on the issue of gap formation, and outline recent progress and 
future points of contact between condensed matter physics and Lattice QCD. 
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1. Introduction 

The recent experimental isolation of single atomic layers of graphite, known as graphene, has 
provided physicists with a novel opportunity to study a system with remarkable electronic and 
many -body properties, which is easy to manipulate experimentally [jl|, §]. Even more recently, 
the advent of experiments utilizing samples of suspended graphene, free from the interference of 
an underlying substrate [Q], has provided unprecedented insight into the intrinsic properties of 
graphene. Among other remarkable discoveries, suspended graphene has been found to exhibit 
a carrier mobility which exceeds that of silicon by an order of magnitude, a fractional quantum 
Hall effect which is indicative of strong electron-electron interactions, as well as a markedly non- 
metallic behavior of the DC conductivity at low temperatures. 

A central property of graphene is that the low-energy electronic spectrum can be described 
in terms of two flavors of massless, four-component fermionic quasiparticles with linear disper- 
sion [Q]. Indeed, due to the hexagonal honeycomb arrangement of the carbon atoms in the graphene 
lattice, the band structure of graphene exhibits two inequivalent (but degenerate) "Dirac cones" 
where the conduction and valence bands touch. Since the energy-momentum relation around a 
Dirac point is linear as in relativistic theories, the low-energy description of graphene bears a cer- 
tain resemblance to massless Quantum Electrodynamics (QED). Nevertheless, an important differ- 
ence is that the Fermi velocity of the quasiparticles in graphene is as low as v F ~ c/300, whereby 
the electromagnetic interaction is rendered essentially instantaneous. However, it should be pointed 
out that such a value of v F is actually unusually large from a condensed matter point of view. 

Such an unusual band structure (often referred to as semimetallic) accounts fairly well for 
the observed properties of graphene sheets deposited on a dielectric substrate. While suspended 
graphene has recently come under intense experimental investigation [gj], its spectrum is yet to 
be computed in a controlled fashion. From a theoretical perspective, the challenging feature of 
suspended graphene lies in the smallness of the dieletric constant e = £ which, in conjunction 
with the small value of v F , results in a graphene analogue of the fine-structure constant a g > 1. At 
such strong coupling, a dynamical transition into a phase fundamentally different from the weakly- 
coupled semimetallic phase of graphene is a strong possibility. In graphene sheets deposited on a 
substrate, such a transition is effectively inhibited due to the suppression of a by the dielectric. 

2. Low-energy effective theory 

The electronic band structure of graphene close to the Fermi level forms the basis of the low- 
energy effective theory of graphene. This band structure is a reflection of the hexagonal arrange- 
ment of the carbon atoms, which can be decomposed into two triangular sublattices A and B. This 
leads to the tight-binding model 

H = -t £ (al i b aJ + h.c.)-t' £ (al 4 a aJ + bl 4 b a j + h.c), (2.1) 

where the operators i {a a ( ) and b^ i {b a ( ) create (annihilate) an electron of spin o at location i on 
sublattices A and B, respectively. The first term (involving t) takes into account nearest-neighbor 
interactions (as first done by Wallace in Ref. |Q]), and the second term (involving t') the next- 
to-nearest neighbor ones. The hopping parameters that give an optimal fit to the experimentally 
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determined band structure of graphene are t ~ 2.8 eV and t 1 ~ 0.1 eV The tight-binding model 
with nearest-neighbor hopping can be generalized to a hexagonal Hubbard model by addition of 
an on-site Coulomb repulsion term, which has been recently studied [Q] within a Quantum Monte 
Carlo approach. 

We shall follow a somewhat different route based on an Effective Field Theory (EFT) treatment 
of graphene [|| g], which has the advantage of describing the physics of graphene directly in terms 
of the relevant low-energy degrees of freedom, namely charged massless fermionic quasiparticles. 
The EFT description of graphene has an additional advantage as it allows for the direct study of 
effects due to the unscreened, long-range Coulomb interactions between the quasiparticles. In 
what follows, we shall formulate a continuum Lagrangian field theory valid at low momenta, much 
smaller than the inverse of the interatomic distance ~ 1 .42 A. 

2.1 Continuum formulation 

The low-energy EFT of graphene may be derived from a tight-binding or Hubbard model 
description augmented by a long-range Coulomb interaction yielding a theory of Nf Dirac 
flavors interacting via an instantaneous Coulomb interaction. The action (in Euclidean spacetime) 
of this theory is 

Nf f If 
S E = - y L] d 2 xdt %D[A ] Va+jglJ dxdt (^ A o) 2 ' < 2 - 2 ) 

where Nf = 2 for graphene monolayers, g 2 = e 2 /e for graphene in vacuum (suspended graphene), 
\jf a is a four-component Dirac field in (2+1) dimensions, A is a Coulomb field in (3+1) dimensions, 
and 

D[A ) = Yo(d + iA )+v F r i d i , i=l } 2 (2.3) 

where the Dirac matrices y, satisfy the Euclidean Clifford algebra {y^,y v } = 28^ v . The four- 
component spinor structure accounts for quasiparticle excitations of sublattices A and B around 
the two Dirac points in the band structure [Q, ^J. The two Dirac points are identified with the 
two inequivalent representations (with opposite parity) of the Dirac matrices in (2+1) dimensions. 
In graphene monolayers, Nf = 2 owing to electronic spin, while Nf = 4 is related to the case of 
two decoupled graphene layers, interacting solely via the Coulomb interaction. Consideration of 
arbitrary Nf is also useful, given that an analytic treatment JTfj| ] is possible in the limit Nf — > oo. 

The strength of the Coulomb interaction is controlled by a g = e 2 /(4nv F e ), which is the 
graphene analogue of the fine-structure constant a ~ 1/137 of QED. It is straightforward to show 
that a g is the only parameter, by rescaling according to 

t' = v F t, A' Q =A Q /v F . (2.4) 



The action (2.2) is invariant under spatially uniform gauge transformations (see Sec. 3.1). Notice 



that since the A field is (3+l)-dimensional, a four-fermion Coulomb interaction of the form 

%(*)%%(*) %&)r Yb&) 



(2.5) 



|x — x | 

is recovered by integrating out A . Nevertheless, for our purposes the original form of the action 



(quadratic in the fermions) as given in Eq. (2.2) is preferable. 
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A central property of the low-energy EFT is that Eq. (2.2) respects a global U(2A^) chiral 
symmetry under the transformations 

Ya -»• exptTTjO,-) Y a (2-6) 

where the matrices IT, are the (2,/Vy ) 2 hermitian generators of \J(2N^), such that for the case of 
graphene monolayers, the group is U(4). It should be noted that the choice of any particular 
representation for the r • is completely arbitrary and is not necessary for any calculational purpose, 
as all relevant information is provided by the Clifford algebra. However, the identification of the 
spinor degrees of freedom with any particular Dirac point and graphene sublattice is dependent on 
the chosen representation. This U(4) chiral symmetry, which strictly speaking is a flavor symmetry, 
can be spontaneously broken down to U(2)xU(2), in which case the excitonic condensate 
acquires a non-vanishing value, signaling the formation of quasiparticle-hole bound states. The 
same group structure is obtained by adding to Eq. (12) a parity invariant (Dirac) mass term 

d 2 xdtm {) %y a , (2.7) 



which breaks the symmetry explicitly. For the extended theory with N* flavors, the symmetry 
breaking pattern is U(2A/y) — >\J(Nf) x\J(Nf). Other symmetry breaking patterns, involving the 
possibilites of magnetic as well as Cooper-like pairing instabilities, have been investigated in 
Refs. [g[llp. 

2.2 Effective action and probability measure 

The partition function corresponding to Eq. ( ^2| ) is given by 

Z = J ®A ®y®y exp(-S E [%, W a Ao\), (2.8) 

where it is possible to integrate out the fermionic degrees of freedom, as S E is quadratic in the \j/ a . 
We thus obtain 

Z = J®A exp(-5|[A ]) det(D[A ]) /v / , (2.9) 

where 

S§ = -L / cPxdt (d,A ) 2 (2.10) 
2g l J 

is the pure gauge part of the action. It is of central importance for the convergence of the Monte 
Carlo algorithm that the above determinant has a definite sign, independently of any particular 
configuration of the gauge field A . One way to establish this property is to proceed by writing 
D[A ] in the form 

d\a]-( m[a ° ] ° )-( M ^ ^ (2 in 

i0] -{ -M[A ])-[ Mt[A ]J' (2 - H) 

where 

M[A ] = cr (d + iA Q ) + v F Oid h i = 1,2, (2. 12) 

which entails a specific choice of Dirac y-matrices. Furthermore, we note that A is real, and that 
the Pauli matrices and the momentum operator are hermitian. The latter implies = —d^, and 
therefore 

det(D) = det(M)det(M t ) = |det(M)| 2 > 0. (2.13) 
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While this property is not affected by the introduction of a parity invariant mass term such as 
Eq. (1.1), the positivity of det(D) breaks down in the presence of a chemical potential. 

The fact that det(D) is positive definite allows for the definition of an effective action that is 
purely real, given by 

S eS [A ] = -N f lndet{D[A ])+S s E [A ], (2.14) 
so that the partition function becomes 

Z = |fM exp(-S eff [A ]), (2.15) 

where P[A ] = exp(— S eff [A ]) > can be interpreted as a positive definite probability measure for 
a Monte Carlo calculation, as outlined in Section |5[ 

2.3 Operator expectation values 

The expectation value of a given operator 0[\(f, ty] dependent on the fermion fields can be 
calculated by taking functional derivatives of the generating functional 

Z[fj,T}] = y n ^4 ^^exp(-5 £ [A ,^,VA,f],T]]), (2.16) 

where source terms have been added to the original action according to 

S E [A ,\j/,y,fi,ri]=S E [A ,\jf,\i/] + J d 2 xdt(\frq+h.c), (2.17) 

such that the modified effective gauge action is a functional of A as well as of the sources T\,f\. It 
is again possible to integrate out the fermionic degrees of freedom and take functional derivatives 
with respect to the sources in the resulting expression 

Z[fj,T]]oc J ^A exp(-5 eff [A ])expf-| d 2 xdtf\lT x \A^rA , (2.18) 

which makes it possible to obtain expectation values in terms of a path integral over A only. While 
this procedure is completely general, it is possible to employ a slightly different approach in order 
to facilitate the computation of the chiral condensate and susceptibility. 

The chiral condensate a, which is the order parameter of the semimetal-insulator phase tran- 
sition in graphene, is defined by 

a = (%Yb): ( 2 -19) 

where the fermion fields are evaluated at the same space-time point. It is useful to note that the 
mass m plays the role of a source, coupled to Y h y h - The expectation value of this operator can 
therefore be obtained by first differentiating the partition function with respect to m and dividing 
by the volume, giving 

° = VzI J d 2 xdt \ff b (x,t)Yb(x,t) exp(-S £ ) = (2-20) 

where a has been defined as an average over the lattice volume occupied by the fermions. On the 
other hand, once the fermions have been integrated out, the derivative with respect to m yields 

° = Vzl ^oT^'N) exp(-5eff[A ]) = i (Tr^ 1 [A ])> , (2.21) 
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where the identities 

det(D[A]) = exp(Tr(log(D[A])), 3det ™) = det(D[A])Tr (V 1 ^] ^) , (2.22) 

have been used. The chiral susceptibility %i mav be found by taking one more derivative with 
respect to m , giving 

da 1 



Xi 



dm Q V - 



Tr 2 (D- 1 )) - (Tr(D- 2 )) - (Tr^ 1 )) 



(2.23) 



which is expected to diverge at a second-order phase transition, and may also yield constraining 
information on the universal critical exponents of the transition. 

3. Graphene on the lattice 



In this section we formulate the lattice version of Eq. Q following Refs. [jD), g]. It should 



be noted in this context that a closely related lattice model of the strong-coupling limit of graphene 
has been considered in Ref. JT4]]. We begin by discretizing the pure gauge sector, where the require- 
ment of gauge invariance implies the use of "link variables" to represent the gauge degrees of free- 
dom. The "staggered" discretization of the fermionic sector is then outlined, as it is the preferred 
choice to represent fermions with chiral symmetry at finite lattice spacing [OL [16|]. Throughout 
this paper, the lattice spacing is set equal to unity, and thus all dimensionful quantities should be 
regarded as expressed in units of the lattice spacing. 

3.1 Gauge invariance and link variables 



The pure gauge part of the Euclidean action, Eq. ( |2.10| ), can be thought of as the non-relativistic 
: of the Lorentz-invariant form \F„ V F^ V where F„ v = d^A v — d v A^, such that 

V MV = F^ + F^ + F^F' = 2F Qj F°J = 2(djA Q ) 2 , (3.1) 



where we have used F^ = (no magnetic field) and d {) Aj = (no electric field induction by a mag- 
netic field), valid in the non-relativistic limit (v F <C c). Thus, for graphene the only non-vanishing 
contribution is the electric field Ej = —djA , which represents the instantaneous Coulomb interac- 
tion between the quasiparticles. 



The gauge action ( |2. 10[ ) is invariant under the time-dependent, spatially uniform gauge trans- 
formations 

A ^A + a(t), y-t-exp^ij dt'a(t')\\i/, (3.2) 

where a(t) is a function of time only. Thus, in spite of its apparent simplicity, the effective theory of 
graphene possesses a truly local gauge invariance, which should be respected by the lattice action. 
To this end, one introduces temporal link variables 

U 0!n = U n = e x V (ie n ), (3.3) 

where 6 n is the dimensionless lattice gauge field evaluated at the lattice point n = (n0)«ij"2>«3)- 
The spatial link variables 

U u = 1 (3.4) 
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are set to unity. It is convenient to express the discretized version of Eq. (2.10) in terms of "plaque- 
tte" variables, denned by 

(3.5) 



^/iv,n ^,n^v,n+e M ^i,n+e v ^v,m 



where, in the present case of a pure Coulomb interaction, the only non-trivial components are U 0! 
and U m . Those plaquette components then correspond to the discretized formulation of the electric 
field. The remaining components corresponding to the magnetic field are equal to unity. These 
statements can be summarized in the expression 

U Liv,n = S^Syi U n U^ +ei + 8 v0 d^ U*U n+ej + 5 m0 5 v0 + S^S vj . (3.6) 

In terms of the gauge link variables and plaquettes, the discretized gauge action corresponding 



toEq. (2.10) is given by [0 



n ji>v 



(3.7) 



where j6 = 1 /g 2 , such that j3 — > v F /g 2 when the rescaling of Eq. ( |2.4| ) is applied. In Eq. (j3~7t), the 
only non-vanishing contributions arise from the terms with (ju, v) = (1,0); (2,0); (3,0); (2, 1); (3, 1) 



and (3,2). Eq. ( p77| ) may be simplified to 



PL 



3 

I 

/=1 



3-£9i (V n ul^ 



(3.8) 



where 9l(x) denotes the real part of x. Eq. (3.8) is referred to as the compact formulation, which 
has been found to pose problems related to spurious monopole condensation in QED and related 
theories [|l^]. On the other hand, the non-compact formulation, which is obtained from Eq. (3.8) 
by expanding ^(U a Ul +e ) to second order in 6, 



giving 



E.N 



1 



1 



+ 



n i=l 



(3.9) 



(3.10) 



is free from such problems [jl^, 19] and allows for a realistic continuum limit. 

3.2 Staggered fermions 

While the discretization of the gauge sector is relatively straightforward, the inclusion of dy- 
namical fermions on the lattice is a notoriously difficult problem. One of the main issues when 
simulating fermions on the lattice is the so-called doubling problem (for an overview, see Ref. Jl7|], 
Chapter 4). This problem is related to the chiral invariance of the fermionic sector, and arises due to 
the appearance of multiple (unwanted) zeros in the inverse propagator. In other words, one is simu- 
lating more fermion flavors than expected, the exact number being dependent on the dimensionality 
of the theory. There exists a number of ways to avoid the doubling problem, but all of them break 
chiral invariance in one way or the other, a fact encoded in the Nielsen-Ninomiya theorem [pC|]. 
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The solution we have chosen for our simulations of graphene is the "staggered" discretization of 
Ref. pT|l. This choice is advantageous for the study of spontaneous chiral symmetry breaking in 
graphene, as it yields the correct number of degrees of freedom while (partially) preserving chiral 
symmetry. The major drawback of staggered fermions is that the full chiral symmetry is restored 
only in the continuum limit, a fact referred to as "taste symmetry breaking". 



In order to discretize the fermionic sector of Eq. (2.2) in a way amenable to computer simu 



lations, a number of choices need to be made. As a first step, fermions are integrated out, and the 



partition function is written purely in terms of the gauge field, Eq. ( |2.15[ ). The fermions are then 
represented exclusively through det(D). One can then attempt to compute the determinant exactly 
for a given 6 configuration, which is feasible due to the low dimensionality of the problem, Alter- 
natively, one may rewrite det(D) in terms of a path integral over complex scalar fields referred to 
as pseudofermions, as is common in Lattice QCD. 



It has been shown in Ref. [ |22| ] that for each staggered flavor one recovers, in the continuum 
limit, two four-component Dirac flavors. Thus, by retaining one staggered flavor, it is possible 
to have exactly eight continuum fermionic degrees of freedom, which is the correct number for 
graphene. The action of a single staggered flavor is given by 

s f E [x,x,e] = - 1 £xnK n , m [e]x m , (3.ii) 

n m 

where the staggered Dirac operator is 

K n , m [e] = l(8 B+eo>m U n - ^m^ + ylli^,"- Ve iim )+™(A, m) (3.12) 



where the phase factors r\ arise from the spin-diagonalization of the Dirac matrices [23]. The op- 
erator K thus replaces D in all expressions for the probability, chiral condensate and susceptibility 
that were derived in the previous sections. As expected from the Nielsen-Ninomiya theorem, the 
staggered lattice action does not retain the full U(4) chiral symmetry of the original graphene ac- 



tion at finite lattice spacing. As shown in Ref. p2], only a subgroup U(l)xU(l) remains upon 



discretization. Spontaneous condensation of %%, or equivalently the introduction of a parity invari- 
ant mass term, reduces this symmetry to U(l). 

Finally, it should be pointed out that the situation concerning graphene is unusually favorable, 
in the sense that the staggered formalism fortuitously provides the correct number of fermionic 
degrees of freedom, as Ny = 2 for graphene monolayers. In general, staggered fermions provide 
only a compromise solution in the sense that some degree of chiral symmetry is preserved, at the 
price of retaining some of the doubling originally present in the discretized fermion action. Indeed, 
if the case of N* = 1 were to be simulated, it would be necessary to resort to the controversial 
"rooting" trick whereby the desired number of continuum flavors is restored by taking the 
appropriate root of the Dirac operator. 

4. Results for j8 c 

The results presented in this work for the non-compact gauge action are partly based on the 



reanalysis of extant data from Ref. [12] and partly on new data on larger lattices. These datasets 
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Figure 1: Result of a simultaneous fit to a (left panel), %, (middle panel) and R (right panel). Red symbols 
indicate data for N x = 28, the remaining ones are for N x — 32. The datapoints for R are grouped according 
to j3 ~ 1 . The fit range has been restricted such that data with m > 0.005 and j3 < 0. 1 are excluded, as data at 
large m and small /3 have substantial finite volume effects. The optimal parameters are /3 6 . = 0.0774(2) and 
8 = 2.93(2), with b = 1.0 fixed. The data satisfy b = 1 very accurately. The error bars on the datapoints are 
obtained via the Jackknife method [E5D. 



correspond to N t =N X = 28 and 32, with the extent of the bulk dimension N b set at either 12 or 32. 
We find that finite-size effects are most pronounced as functions of N t and N x at small /3, while the 
effect of N b on the results is negligible. Our analysis proceeds by first determining the condensate 
a, susceptibility %i an d R = Xi m o/ a as a function of j8 and m . The second step of the analysis 
consists of fitting an equation of state (EOS) to a and Xi m order to obtain estimates of the critical 
coupling j3 c . and the critical exponents. Our EOS is of the form 

mo X(P) = Y(l3)o b + o s , (4.1) 

where the exponent b = jS — 1/5. This form has been used previously in the context of Lattice 
QED in (3+1) dimensions [|2(]]. Here, we have referred to the critical exponent using the notation j8 
to avoid confusion with the inverse coupling j3 . The dependence on the critical coupling f5 c enters 
through the expansions 

£)+..., 

where terms up to X l and Y { have been retained. Higher-order terms were found to have a very 
small effect on the analysis, and have thus been discarded. 

We find that a simultaneous fit to all data for a, Xi an d R yields the most stable results. It 
is also necessary to carefully consider finite-volume effects and the impact they have on the fit 
parameters. We have therefore excluded datapoints with ra > 0.005 and /3 < 0.1. In this way, we 
find consistency with mean-field exponents, 5 = 3 and 6 = 1. Including data at lower /3 and larger 
m suggests 8 ~ 2.2, however such fits have a much higher x 2 a nd compare unfavorably with the 
data on R. We find that N x = 28 and N x = 32 give consistent results using the restricted dataset. 



X(fi)=X Q +X l [l 



+ ..., 



(4.2) 
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5. Experimental situation 



We now turn to the question whether experiments which measure the DC conductivity of sus- 
pended graphene provide any evidence for semiconducting behavior which would follow naturally 
from the excitonic gap scenario. While a full LMC calculation of the conductivity is not yet avail- 
able, a simplified analysis in terms of a Kubo description of gapped quasiparticles has recently 
been given in Ref. [27], where the data of Ref. [ |2"8| ] on the suspended graphene devices "SI", "S2" 



and "S3" were analyzed in terms of the expression o = o q + o bg . Here o q is the quasiparticle 
contribution intrinsic to graphene, while the "background" component a bg is device-dependent. 

The Hamiltonian describing Dirac quasiparticles with a gap A and Fermi velocity v F ~ c/300 
is given by H = a l v F k l + o 2 v F k 2 + o 3 A/2, where the a, are Pauli matrices. The contribution o q of 
the Dirac quasiparticles to the DC conductivity of a graphene monolayer is then 

„ t = gf r de r« WM /P«-¥-<w-/«.« + s-fto (5 „ 

' n 2 J-oo J a/2 CO 

where /3 = \/k B T , the Fermi function is given by f(x) = 1/(1 +exp(x)), pL denotes the chemical 
potential and the factor of 4 accounts for the spin and valley degrees of freedom. Then 



+ 



<g 2 + A 2 /4 
£ 2 -A 2 /4 



CO 



CO 



CO 



CO 



CO 



CO 



















(5.2) 



where tj is the scattering rate of the quasiparticles, which can be accounted for Q29Q by broadening 
the delta functions according to nd^ (x) = f] / '(x 2 + T] 2 ) . In the DC limit, the integral over £ in 
Eq. yields 

r°° 1 A 2 — 4l7l 2 

dS ^ Q ^,e) = — - - r7 -J^arg (A 2 -4z 2 ) , (5.3) 



/A/2 " " ~'~ 2X 167T£T] 

where z, = £ + ir\. 

The inclusion of the background component a bg is motivated by the fact that the minimal 
conductivity in graphene samples is non-universal. In suspended graphene, it is much smaller than 
in graphene samples on a substrate and furthermore strongly sample-dependent. There is also a 



clear tendency of the minimal conductivity to decrease with increasing sample purity jpOQ. It is 
likely that the minimal conductivity in graphene is formed of several components, including the 
effects of impurities ^T]], inhomogeneity [|32|] and effects due to invasive metallic contacts [33]. In 
our analysis, we have used the phenomenological form 



o bg = o exp[-(T /T) a ], 



(5.4) 



which allows for the slight empirical T-dependence of a, . The empirical data of Ref. [El] is shown 



bg- 

in Fig. ||, together with fits in terms of o q and o hg . It is noteworthy that the data display a distinct 
"knee" at ~ 30 K, which in terms of the present description is interpreted as the temperature below 
which thermal activation is negligible. Thus, in order to determine a bg in an unbiased fashion, we 
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Figure 2: Left panel: quasiparticle (QP) and background (BG) components of a{n = 0, T) for the suspended 
graphene devices S1-S3, as determined in Ref. [p7[|. The empirical data is reproduced from Ref. [pq]. 
All devices show a "knee" separating thermally activated and background regions. Right panel: a(n,T) 
determined from a fit to o(n = 0, T) and the resistivity p (n, T), reproduced from Ref. p^. 



first fix a and T using data in the extreme low-J region. One may then subtract o hg at all T , 
and determine fir] and A by fitting o q to the resulting dataset. A simultaneous fit of o q and o bg 
confirms the validity of this procedure. While v\ (T) is a priori unknown, a scenario of constant /3tj 
is strongly favored by the available data in the range 35K<T<150K. 

Our findings in Ref. [|7|] suggest that the suspended graphene devices of Ref. [ 28 ] exhibit a 
thermally activated conductivity a q , which is well described by Eq. ( |Q| ) from T ~ 150 K down to 
T ~ 35 K, where the signal is lost due to limited measurement accuracy. The determined bandgaps 
are in the range A ~ 25 —40 meV, whereas all samples were found to favor j3rj ~ 0. 1 indicat- 
ing a scattering rate which increases linearly with T. A natural scattering mechanism with such 
properties is provided by the long-range Coulomb interaction [34] up to logarithmic corrections. 
Furthermore, this value of j8 f] is consistent with the high carrier mobilities and long mean free 
paths reported in Ref. [e3], as well as with theoretical estimates of the mean free path due to long- 



range Coulomb scattering [34]. Specifically, for T = 35 — 150 K we find r\ = 3.5 — 15 K, with 
corresponding mean free paths of Tiv F /r\ ~ 0.5 — 2.0 jUm. Fits with no gap (A = 0), constant T7, or 
zero background were found to be incompatible with data. 

As shown in Fig. |2|, these conclusions are consistent with the observed o(n), which depends 
sensitively on the value of /3tj determined from data at n = 0. Furthermore, the interpretation 
of the observed c(n,T) as due to thermal activation accounts, in a natural way, for the observed 
transitional density n* above which o{T) reverts from insulating to metallic. The determined 
scattering rate f]{T) is also suggestive of long-range Coulomb scattering, which is consistent with 
the ultrapure character of the suspended graphene samples of Ref. [28]. 
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6. Conclusions 



We have reviewed the Lattice Gauge Theory approach to the low-energy EFT of graphene, 
with the aim of introducing this technique to a wider audience and motivate the application of 
this approach to systems beyond monolayer graphene. Our calculations within this the graphene 
EFT indicate that it displays a chiral phase transition at a critical coupling of j8 c = 0.0774(2), with 
critical exponents that appear consistent with mean-field theory. Spontaneous chiral symmetry 
breaking in the graphene EFT would lead to the appearance of a gap in the quasiparticle spectrum, 
directly linked to the formation of quasiparticle-hole pairs (excitons). It is conceivable that such a 
transition occurs in suspended graphene, where the strength of the Coulomb interaction attains its 
maximum value. 

In an effort to clarify whether currently available experimental data on suspended graphene 
provide any evidence for the excitonic scenario, we have reviewed the status of such measure- 
ments, which show a definite (though relatively mild) insulating trend at low T in the vicinity 
of the neutral point. We have presented an interpretation of the observed anomalous temperature 
dependence in terms of the excitonic gap scenario, and tentatively found that the data may be 
consistently explained in terms of gapped Dirac quasiparticles (A ~ 30 meV) with the long-range 
Coulomb interaction as a natural candidate for the dominant scattering mechanism. Further experi- 
mental studies of the conductivity at low T in suspended graphene are clearly called for, preferably 
minimizing the effects of invasive metallic contacts. 

Further investigations using the Lattice Gauge Theory approach are in progress, including 
the renormalization of v F due to the Coulomb interaction, the magnetic catalysis of a semimetal- 
insulator transition at large external magnetic fields [35, 36], and the critical temperature for exci- 
ton condensation in graphene bilayers [37]. The computation of transport properties involves the 
extraction of spectral functions in Minkowski spacetime from Euclidean time lattice data. Such 
calculations are also feasible nowadays, as Bayesian analysis coupled with the Maximum Entropy 
Method has been successfully applied to QCD []3"8l]. Areas of interest include the electrical con- 
ductivity and viscosity | |39| ] of the quasiparticles in graphene. Due to the flexibility of the Lattice 
Monte Carlo and Quantum Monte Carlo methods, further applications are likely to arise in the 
closely related fields of graphene nanowires, high-r c superconductors and hexagonal optical lat- 
tices, to name a few. In summary, the application of Lattice Gauge Theory to condensed matter 
problems appears poised to develop into a highly fruitful field of study. 
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